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It has recently been shown that F-theory based constructions provide a potentially promising avenue 
for engineering GUT models which descend to the MSSM. In this note we show that in the presence 
of background fluxes, these models automatically achieve hierarchical Yukawa matrices in the quark 
and lepton sectors. At leading order, the existence of a U(l) symmetry which is related to phase 
rotations of the internal holomorphic coordinates at the brane intersection point leads to rank one 
Yukawa matrices. Subleading corrections to the internal wave functions from variations in the 
background fluxes generate small violations of this (7(1), leading to hierarchical Yukawa structures 
reminiscent of the Froggatt-Nielsen mechanism. The expansion parameter for this perturbation is 
in terms of y/acuT- Moreover, we naturally obtain a hierarchical CKM matrix with V\2 ~ V21 ~ e, 
0^ ' ^23 ~ V32 ~ e^, V13 ~ V'ii ~ where e ~ y/otcuT, in excellent agreement with observation. 
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Flavor physics remains a poorly understood aspect of string based constructions as well as phenomenological models 
which aim to reproduce the Standard Model or MSSM at low energies. In this note we show that the minimal form of 
the F-theory GUT models recently developed in [l|,[3,[3 (see also 0, d, d, 0, B i, [iSl ) automatically generates viable 
Yukawa matrices in the quark and lepton sectors. In the present context, a minimal implementation of an F-theory 
T ' GUT simply reflects the discrete choice in the geometry that we consider models with the fewest possible number of 
matter curves and intersection points necessary for compatibility with the MSSM. 
The MSSM superpotential contains the terms: 



Wmssm = ■ Q'Um^ + X'J ■ QW^d + XI' ■ L'E^Ha + ■■■ (1) 

m 

' where i and j are generation indices so that refers to the right-handed top quark, and the A's denote 3x3 matrices. 
Although genericity arguments might suggest all of the associated masses should be comparable, there is a well-known 
^ hierarchy 

^ (m„,TOc,mO (0.003, 1.3, 170) X GeV (2) 

[md, ms,mb) ~ (0.005, 0.1, 4) x GeV (3) 
(me, m^, TO^) ~ (0.0005,0.1,1.8) X GeV. (4) 

O 

To leading order, the inter-generational ratios of masses between particles with the same gauge quantum numbers are 
insensitive to running effects. On the other hand, running to the GUT scale reduces the heaviest generation quark 
masses by roughly a factor of three. 
^ Letting denote unitary matrices such that V^XuV^'' and Vj^XdVf ^ are diagonal, the norm of elements in the 

CKM matrix Vckm = Vu^d^ are also hierarchical [llj |: 

0.004 

\VcKM{M^eak)\ ^ | 0.23 0.97 0.04 I . (5) 

0.99 

This is encapsulated in the Wolfenstein parameterization of Vckm UM ■ 

A AX^ {p - iT]) 

Vckm = | -A 1 - AV2 AX^ | + OiX") (6) 

^AA2 1 
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where A ~ 0.23, A ^ 0.82, p ~ 0.22, i] ~ 0.34 at the weak scale. To one loop order, only A evolves with scale [Tsl. 
In fact, even including two loop running effects, A only changes by an order one factor [TH - 

This type of hierarchy is potentially present in models with a global U{1) flavor symmetry as in the Froggatt- 
Nielsen mechanism [l^. For example, assigning FN charges 4,2,0 to the three generations of Q, U superfields and 

introducing a chiral superfield Xp^ with FN charge —1, the superpotential term (Xfat/Mfat)^^'^"'''^^^'^""'' • Q^U^ 
will induce hierarchical Yukawa couplings once Xp^ develops a suitable vev. Similar considerations apply for the 
remaining Yukawa matrices. 

On the other hand, the three generations all possess the same gauge quantum numbers with respect to the Standard 
Model gauge group. For this reason, introducing an horizontal symmetry may appear somewhat ad hoc. In this note 
we show that an approximate symmetry of the local geometry plays a similar role to that present in the Froggatt- 
Nielsen mechanism. Nevertheless, we also find that the exact structure of the Yukawa matrices differs from single 
field Froggatt-Nielsen models due to the profiles of the internal wave functions. 



II. PARTIALLY TWISTED GAUGE THEORY AND F-THEORY GUTS 

In F-theory GUT models (see [UllBl for reviews), the gauge fields of the MSSM descend from the cight-dimcnsional 
worldvolume of a seven-brane wrapping a complex del Pezzo surface S with a GUT gauge group. This GUT group is 
broken to SU{i)c x SU{2)l x U{1)y via an internal flux through the seven-brane in the U{1)y direction of the GUT 
group [3, Q. The chiral matter of the MSSM descends from zero modes of six-dimensional fields localized on matter 
curves in S. In the classical limit, the zero mode wave functions vanish off of the matter curve. Strictly speaking, this 
is only approximately true because these wave functions have non-vanishing support in directions transverse to the 
matter curve. In purely gauge theoretic terms, these wave functions can be derived from an appropriate internal local 
Higgsing of a parent eight-dimensional theory [205. this parent theory, there is an adjoint- valued field (j)- When <j) 
develops a vev which breaks the parent theory gauge group, adjoint- valued fermions in the parent eight-dimensional 
gauge multiplet descend to "bifundamentals" which are trapped along complex codimension one matter curves defined 
by the vanishing locus = 0, much as in the Nielsen-Olesen vortex jl9|. The trapped zero modes obey the system of 
equations ^li]: 

UJS A dA^Ja + 2 [01 Xa] 
dAXa - [0, V'a] 

where tus denotes the Kahlcr form on S, and Xa as well as ipa denote fermions in the 8d gauge multiplet. The subscript 
A reflects the background internal gauge field from the GUT seven-brane, as well as potentially other seven-branes 

which intersect the GUT seven-brane. In addition the "O ^ ^^^^'^ ^ " reflect the possibility of contributions from higher 

dimension operators induced by other background fluxes of the compactification associated with more general p-form 
potentials. Here, M^, is the characteristic mass scale of the F-theory compactification which is related to the GUT 
scale Mqjjt and acuT as 01 : 

= ao'uT ■ MhuT- (9) 

In general, the effects of these background fluxes can non-trivially mix with such higher-form potentials. For example, 
in an abelian gauge theory on a D-brane, the presence of an NS B-field shifts the field strength as: 

F'=F + B. (10) 

For simplicity of presentation, consider the case where the background field strength is parallel to the matter curve. 
Letting z±^ denote the coordinate normal to the curve = 0, and z the coordinate on the curve = 0, the trapped 
zero mode wave function is of the form: 

* - /(z, z) exp(-7z_LzT), (11) 

where / is a 9a zero mode of the bundle defined by the gauge field along the curve = 0, and 7 depends on the 
overall scale specified by the Kahler form. Note that the ratio of any two zero modes of Oa can be expressed purely 
in terms of holomorphic functions of z, a fact that will be crucial in the analysis to follow. Further note that by 
rescaling the Kahler form so that 7 ^ 00, the choice of representative used in the topological field theory can be made 
arbitrarily peaked along the matter curve. 



+ (7) 



+ f^l) (8) 
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The intersection of three matter curves in S leads to Yukawa coupHngs among the corresponding matter fields. The 
Yukawa couplings between the Higgs fields and chiral matter are given by triple overlap integrals of the form: 

= J A**$^'. (12) 
s 

Here, A^'*$-' may be viewed as shorthand for the triple product of the wave functions A/^^^g^E'^, A/^^^'q^-^ or 

A/z^^'^^E'^ in the obvious notation. In the limit where the internal field variations (other than the Higgs) is constant, 
this integral reduces to a product of the three wave functions at the various intersection points of the matter curves: 

X^^ ^J2^ip)^\p)'f\p) (13) 

p 

where p denotes a point of triple intersection. Indeed, this superpotential is the descendant of the superpotential for 
bulk modes in the parent theory which makes no reference to the Kahler form. As such, we are free to rescale the 
Kahler form dependence until the overlap of the corresponding wave functions is only non-zero at the point of mutual 
intersection. 

This is quite analogous to the computation of perturbative Yukawa couplings in the context of (2, 2) sigma models, 
where in the classical limit the Yukawa couplings are given by classical intersection theory. In the A-model, this result 
is then deformed by worldsheet instanton corrections and in the B-model by the variation of the complex structure. 
The setup in F-theory is parallel to the B-model setup and we can view the above Yukawa coupling as the analog of 
the limit of "large" complex structure. Just as in the B-model at finite complex structure, there will also be subleading 
corrections to equation (|13p . Such corrections can occur as a result of higher dimension operator contributions to the 
effective action, which will deform equations ([7]) and ([5]) . 

From a bottom up perspective, a minimal implementation of the chiral matter of an SU{5) GUT requires one 
matter curve for the 10m and one for the 5m where SU{5) respectively enhances to 5*0(10) and SU{6). The Higgs 
5h and 5h localize on distinct curves where SU{5) enhances to SU{6). The background flux through these curves 
dictates the number of chiral matter zero modes in the low energy theory, while a suitable background flux through 
the Higgs curves achieves doublet-triplet splitting 0- 

The superpotential terms 5h x Wm x 10a/ and 5h x 5a/ x 10m respectively originate from enhancements to Eq and 
50(12). Geometrically, this requires the 10a/ curve to intersect itself, as well as the 5h curve at a single point, and 
for the 10m curve to intersect the 5m and 5h curves at another point of 5 We shall respectively label these triple 
intersection points as Pup and Pdown- While additional intersections are in principle possible, the minimal, and most 
generic situation will typically only contain one intersection point of each type. Note that in principle there can be 
additional points where the singularity type enhances since the matter curve may consist of more than one irreducible 
component. In the special case where the wave functions are given identically by a sharply peaked Gaussian which 
is concentrated along a matter curve, a single triple intersection for each type of Yukawa leads to rank one Yukawa 
matrices [3]. For example, the up quark Yukawas are given by: 

= ^nAPupWQMKiPup) (14) 

which is manifestly a rank one matrix. If there were no further corrections, the resulting theory would contain a 
massive top quark, and massless up and charm quarks. Similar considerations hold for the down type and lepton 
masses. 

In evaluating subleading corrections to the Yukawas, it is more natural to go to a canonical orthonormal basis of 
wave functions where this rank one structure in the Yukawa matrix is more manifest. To this end, let zi and Z2 denote 
two local coordinates on S such that zi = and Z2 = denote the matter curves where the wave functions $^ and 
of equation respectively localize. By an appropriate unitary change of basis, we can always arrange for ~ zj*"' 
and $^ ~ Z2~'' along the respective curves where the wave functions localize (in particular we can choose these wave 
functions to be orthonormal by allowing subleading pieces in the powers of Zi in each wave function). In this basis, 
A'^^ ^ whereas all the other components of A*-' vanish. In terms of the wave function overlap integral of equation 
this vanishing structure in the Yukawa matrix can be ascribed to the presence of a geometric U{l)i x U{1)2 
action on the Zj: 

Zj ^ exp{iaj)zj. (15) 
The terms of the integrand which vanish are precisely those which are not invariant under this [/(l)ixt/(l)2 symmetry 
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III. WAVE FUNCTION DISTORTION AND YUKAWA HIERARCHIES 



We now compute subleading corrections in the Yukawa matrices due to distortions in the profile of the internal 
wave functions. In principle, such distortions from various fluxes induced by p-form potential of the compactification. 
The ones of relevance for distorting the Yukawas correspond to higher dimension operator deformations of equations 
d?]) and ([8]). The overlap integral of interest is given by: 



X'' = Jd'z,d^z,-niz,,zT,Z2,z^)- (j^y (S)' ' ^^^^ 



where i?i and R2 denote the characteristic lengths of the curves zi = and Z2 ~ 0, and 17 (zi,zr, Z2,'Z2) is a function 
which includes the contributions from the Higgs field wave function, as well as other non-holomorphic terms associated 
with the other wave functions. Note that all information about flavor is contained in the term zf~^Z2~'' ■ 

As noted above, a non-zero contribution to A*-' can only occur when a contribution to the integrand is invariant 
under the local U{l)i x C/(l)2 action of the coordinates. A non-zero entry for the A^-' when i and j are both different 
from 3 therefore requires some contribution from to carry non-trivial U{l)i x U{1)2 charge. 

Such corrections originate from general fluxes of the compactification. Since is holomorphic and in the local 
geometry lies in a U{1) x U{1) subgroup of the parent local gauge symmetry, without loss of generality we can choose 
local coordinates so that 0i = zi and (/<2 = Z2 in the neighborhood of the intersection point. Our aim will be to 
characterize possible violations of the U{l)i x U{1)2 symmetry due to background fluxes. Owing to the symmetries 
of equations ([7]) and ([8]) and their deformation by possible higher dimension operators, the leading order distortion 
of the wave functions from background fluxes is a linear combination of the U{l)i x U{1)2 invariants ZiZj so that: 



where ^''-'^^ denotes the contribution to the wave function in the limit where the field strength vanishes, and Ai^^ is a 
2x2 matrix. Although the off-diagonal contributions to this quadratic form indeed violate the "axial" combination 
of J7(l)i and U{1)2, such terms cannot contribute to the Yukawa couplings [2^. 

We stress that the distortion due to M can in principle originate not just from gauge fields, but from more general 
fluxes induced by p-form potentials of the compactification, which will also mix with these gauge fields. Indeed, while 
gauge field flux alone does not turn out to distort the Yukawa couplings, other fluxes do [26]. Further note that by 
dimensional analysis, the overall scaling of Ai is fixed as A4 ^ Mqjjj,. 

More generally, the field strengths will vary over points of the geometry. In an adiabatic approximation we can 
simply retain the form of ^' but with M''^ now given by a function of zi,zr, Z2 and zj. We now perform a series 
expansion in derivatives of around the point of triple intersection: 

M'^z.zj^ 2^ \L , , ■z'^zi'z^Z2'' -z.z,. 18 



The crucial point is that the monomial z^zi'z^zi'"' has U (l)i x 17(1)2 charge {k — l,m — n). Expanding the exponential 
of equation (|17p . it now follows that there will generically be corrections to the rank one Yukawa coupling. 

We now compute the leading order behavior of the Yukawa matrix entry A'^ by performing a series expansion of f2 
in equation (|16p . U{l)i x U{1)2 charge conservation implies that the only terms from a series expansion of CI which 
can contribute are of the form: 

{z,zTr{z2Z^)''-izT)'-'{z^f-'. (19) 

Each power of z reflects the presence of a derivative in the series expansion. 

There are in principal two expansion parameters available which will lead to different hierarchies in the Yukawa 
couplings. The flrst, and perhaps most obvious possibility is that we can simply count the number of total z's which 
appear in a given term. This amounts to performing a derivative expansion in the gauge field strength. In this 
case, the leading order contributions will always originate from terms where a = b = 0. The contribution from the 
derivative expansion is given as: 

m!-^-^ [rimJ [r2mJ ^^^> 
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where A4 is shorthand for the generic entry A^*^ , c*-' is an order one coefBcient which reflects the evaluation of the 
integral in equation (fT6|) . 

On the other hand, there is another natural expansion parameter given by expanding in successive powers of the 
first gradient of the flux. The leading order perturbation to the Yukawas from this class of terms is: 

*A?„^c«^f%^y'7¥lfV". (21, 



To estimate the magnitude of each 5X^^ , we note by dimensional analysis, A4 ~ Mqjjj,. Further, each successive 
derivative of M contributes an extra factor of l/i?i or 1/R2- Setting: 

eo = M{0)/Ml ei = {M^Ri)-^ , £2 = (M,i?2)"' , (22) 

we thus obtain: 

5\%EB.^^'' -'^o-isif-'is^f-' (23) 

54lx ^ ■ i^^^if^' {^oe2f'' (24) 

adopting a parametrization where £1 ^ £2 ^ n and ^ eo£2, the Yukawa coupling matrix therefore assumes the 
form: 

/ \ / ^3 • £2 k2 . g2 ^ . g2 \ / £8 ^6 g4 

A = A(o) + SXder + SXflx - + • ^ • £2 ^2 + ^6 ^4 ^2 . (35) 

V 1 / \k-£^ £2 J \£^ £' 1 

where each entry is multiplied by an order one constant. At this point, it is important to note that although the 
analysis crucially depends on the violation of a U{1) symmetry, the structure of SXder is somewhat different from the 
simplest Froggatt-Nielsen parametrization mentioned earlier. Note, however, that SXp^x is indeed consistent with a 
single field Froggatt-Nielsen model. 

In the geometry, there are in fact three k's and three e's because there are three distinct Yukawa matrices of 
interest. Indeed, the background hyperflux will generate order one distortions in the lepton doublet and down quark 
wave functions so that the entries of A/ and Ad can differ by order one constants, so that there is no a priori mass 
relation between the down type quarks and the charged leptons f2:] . 

Because Ri ^ Af^^y, one might at first think that since e ~ k ~ \/0iGUTi that the DER expansion will always 
dominate over the FLX expansion. However, there can be a further enhancement in the FLX expansion because it 
is more strongly on how chiral matter couples to background fluxes. Since ^^ogut is not extremely small, the FLX 
expansion can in principle dominate over the DER expansion in Yukawas involving fields which have large couplings 
to the background flux. 

Determining whether SXdeb. oi" SXpLX dominates at each order depends on how the left and right chiral matter 
wave functions couple to the background fluxes. In F-theory GUTs, a background flux in the hypercharge direction 
is always present [2, 9], so to illustrate this point, consider the coupling of the matter wave functions to this Hux. 
We shall assume that the strength of this coupling to background fluxes holds for other more general fluxes since in 
general, these higher form fluxes can mix with the gauge fields, as in equation (jlOp so we shall use it as a rough guide 
to determine which type of expansion dominates for a given type of particle species. In the case of hyperflux, the 
internal wave functions couple to this particular background flux in strength proportional to their hypercharges. It 
follows that the hyperflux contribution to Ai''^ is proportional to the corresponding hypercharge. The term SX^pj^^ 
can dominate over (5A^^j^ provided: 

Y .^V>y .^.r^V (26) 

where F denotes the generic field strength, and l^ax denotes in an integral normalization the maximal norms of the 
hypercharge for the left and right chiral matter wave functions associated to each Yukawa. In particular, we have: 





lepton 


up 


down 


y 

max 


6 


4 


2 



(27) 



Using the representative scaling F/M^ ~ 1/ {RM^)^ ~ (^\iuT yields: 

.F=r„,ax-aJ/^T^l. (28) 
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1/2 

where T measures the relative strengths of the flux and gradient expansions. Plugging in aj^y 
integral values for the hypercharges yields: 





Icpton 


up 


down 


T 


1.2 ~ 1 


0.8 - 1 


0.4 < 1 



0.2 as well as the 



(29) 



This rough computation suggests that for the down type Yukawas, SXder. dominates, whereas for the up type and 
charged lepton Yukawas, SX'pj^-^ could potentially dominate. Even though this is by itself quite heuristic, we will 
shortly find that comparison with observation indeed corroborates this picture. 
Setting /t ~ e the Yukawa matrices for Xder and Xflx are: 



^DER • 







£6 




, Xflx ^ 


£6 


£4 


£2 






£2 


1 



this leads to distinct mass hierarchies: 



Xder toi : m2 : ms £^ : £^ : 1 
Xflx ^ mi : m2 : ^ : : 1. 



Next consider diagonalization of A via the matrices and . In fact, it is enough to diagonalize AA^: 

, (Xflx) (Xflx)'' 



( £5 £3 



(Xder) (Xder) 



( £6 p4 



£3 £2 1 



(30) 



(31) 
(32) 



(33) 



£4 £2 1 



Note that these matrices have a structure of the form £°-i^°-i which is consistent with a single field Froggatt-Nielsen 
model with charges 3,2,0 for {Xder) (Xder)^ and 4,2,0 for (Xflx) (Xflx)^ ■ In this case, the matrices and 
are of the form gki^'^il so that: 



^DER 



yR 



DER 




1 ^FLX ^FLX 



( 1 



1 £2 
£2 1 



(34) 



In the above we have focussed on the overlap of the wave functions at a single point of the geometry. At each 
intersection point, there is a local choice of f/(l)'s with respect to which the given mass hierarchy of the generations 
is manifest. Indeed, this is a basis independent statement. 

On the other hand, the CKM matrix Vckm = V^V^'^ measures the mismatch between the matrices which di- 
agonalize the Yukawas in the same basis. Preserving the hierarchy of the CKM matrix requires that the hierarchy 
present in and Vj^ is manifest in the same basis for the Q wave functions. For this to be the case, the eigenspace 
decomposition defined by the local C/(l)'s at the two different interaction points must remain relatively unchanged. 
This occurs provided the distance between the intersection points Pup and Pdown (both of which lie on the Q matter 
curve) should be less than 0.1 x MJ^y. Note that this is to be compared with the length scale of the curves, which is 
roughly Mq^j,. While this is certainly possible, and involves only a very mild fine tuning in the location of the two 
points of triple intersection, it is suggestive of a higher unification structure which could naturally accommodate the 
strict identification p^p = Pdown- The presence of a higher unification structure is also in line with the appearance 
of a U{l)pQ factor in F-theory GUT scenarios which incorporate the effects of supersymmetry breaking [Sj]. For 
the purposes of this paper, we shall only assume that pup and Pdown are sufficiently close that this subtlety can be 
neglected. 

Under this mild assumption, it is in fact now possible to determine the form of the CKM matrix. When either V'" 
matrix is determined by the derivative expansion, k ^ e dominates over order £2 terms so that: 



/ 1 



F-th 
CKM 



£ 
1 

■r2 



/ 



1 



1 

1/2 



GUT 
3 /2 

V '^GUT O^GUT 




(35) 



By inspection, the 3, 2, Froggatt-Nielsen charges dominate in the CKM matrix. In addition, the gap in the Froggatt- 
Nielsen charges can be traced back to the tensor structure of the fiux, which distinguishes it from a scalar. In the 
less likely possibility (which is indeed not realized) where both V^^'s are given by Vpj^xj the form of Vckm is again 
of the form Vp^^. 
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IV. COMPARISON WITH EXPERIMENT 



In the previous section we obtained a parametrization of the up, down and Icpton Yukawas in terms of the parameters 
e and n. In this section, we show that the most natural estimates for e and k dictated by the GUT structure are in 
beautiful accord with observation. Using the explicit values: 



a. 



1/2 
GUT 



0.2, 



we can now estimate the CKM matrix to be: 

/ 



F-th 
CKM 



1 



1/2 



1 



1/2 
^GUT 
3 /2 



3/2 
^GUT 

OLGUT 

1 



^ 1 0.2 0.008 
0.2 1 0.04 
^ 0.008 0.04 1 



(36) 



(37) 



Comparing with the measured values [ll| : 



\VcKM{Mweak)\ 



' 0.97 0.23 0.004 
0.23 0.97 0.04 
^ 0.008 0.04 0.99 



(38) 



reveals a beautiful match between theory and observation! Even though we should only have expected an agreement 
in parameters up to order one effects, we find this close match with the data reassuring. 

We also find close agreement with the expected mass hierarchies in the quark and charged lepton sectors. Due to 
the different structure of the perturbations 6\der and 6\flx, there are in principle two different hierarchical mass 
ratios given by: 



^DER ^ mi : 1712 : ^ e'^ : : 1 



X 



FLX 



mi : m2 '■ 



-4 . 



1. 



(39) 
(40) 



As discussed previously, the dominant contribution to is expected to be from SXflx- On the other hand, we 
have observed that the dominant contribution to A„ and A; could in principle be either SXjjer or SXflx- 

To compare with observation, we can fit the observed masses to a given value of e. At the present crude level of 
analysis, it is enough to neglect the effects of running in such a match because we will be taking third and fourth 
roots of order one numbers. It turns out that the best match to our scenario is achieved when both A^ and A/ are of 
the latter type with corresponding values for e: 





IpLX 


uflx 


doER 


e 


0.36 


0.26 


0.27 



(41) 



1/2 

which are all within order one factors of oiq^jj, ~ 0.2! Fixing the values of the top, bottom and tau mass to their 



observed values, we can now use the ratios of ([39)) and ((40)) to extract the masses of the lighter generations. Comparing 
the F-theory result with the observed values yields: 



F-th\ 



[md,m^ 
{me,m. 



F-th 



F-th\ 



) , [ms,m 
m^,,m:^ 



{mb,\ 



F-th\ 



F-th\ 



(0.003, 0.004), (1.3, 0.8), (170, 170) x GeV 
(0.005, 0.006), (0.1, 0.08), (4, 4) x GeV 
(0.0005, 0.0005), (0.1, 0.03), (1.8, 1.8) x GeV, 



(42) 
(43) 
(44) 



which by inspection, are quite similar. 

The parametrization of the up and charged lepton Yukawas in terms of SXder does not reliably reproduce these 
same mass ratios. One way to see this is to restore the explicit e and k dependence in the mass ratios. There are two 
independent ratios of masses, so in this case we can exactly solve for e and k with the result: 





Ider 


uder. 


doER 


£ 


0.88 


0.40 


0.33 


K 


0.07 


0.05 


0.22 



(45) 



1 /2 

Note that the e's for Ider and uder are both greater than Uq^j.^ whereas the k's are both smaller. Indeed, the 
ratio e/ K. for these two cases is now an order ten number, rather than the order one parameter expected from general 
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considerations. Note, however, that the ratio e/k for cIder is indeed an order one number, which is in accord with 
our scenario. 

The match to the CKM matrix and the mass parametrization of table (|4T|) lends considerable credence to the 
simple picture of wave function distortion we have found. We note that although wc have implicitly worked within 
the framework of the MSSM, the hierarchies we have found are more general, and only require a supersymmetric 
structure near the GUT scale. 

Up to now, we have only discussed the hierarchy in the mass ratios. In fact, our model also predicts the actual 
masses for the top, bottom and tau mass as well 



cny 
If: 



m 



^-'''{Mgut) ~ aJuj, ■ (Hu) = aJuT ' (^> sin/? (46) 



mj^-'\MGUT) - m^~'\MGUT) = c?^tjT ' (^<i> = 4aT ' (^> cos/3. (47) 

Using the value (iJ) ^ 170 GeV, and assuming the large value tan/3 ^ 30 common to the models discussed in Q 
yields: 

m[-*'\MGUT) 20 GeV (48) 
mj^-"\MGUT) - m^'^'iMGUT) - 0.6 GeV, (49) 

which are to be compared with the observed values run up to the GUT scale. These running effects are roughly given 
by a factor of three reduction for the quark masses so that: 

mt{MGUT) - 55 GeV (50) 
mb{MGUT) - niriMGur) 1-5 GeV (51) 

which matches the values of our scenario up to order one factors (in fact, a factor of three). In tandem with our 
estimate for the mass ratios, we thus obtain a crude estimate for the masses of all the quarks and charged leptons in 
terms of (H), tan/3 and ugut- 

V. DISCUSSION 

In this note we have shown that simply demanding a local F-theory GUT model with the most generic features 
automatically implies hierarchical structure in the Yukawa couplings. Moreover, we have found that simple estimates 
on the behavior of the wave functions near points of triple intersection are in remarkable accord with observation! 
While it is certainly quite satisfying to see this type of structure naturally emerge from a string based model, it is 
also possible to go further. 

The first point is that the general form of the matrices we have found can in principle be extended to an arbitrary 
number of generations. Indeed, it is interesting to ask whether this model constrains the existence of additional 
generations of quarks [23| . For example, in a four generation model, the analogue of equation Ij35p is now determined 
by the FN charges 4, 3, 2, 0, leading to: 
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(52) 



Mixing in the three lightest generations corresponds to the upper left 3x3 block of this matrix, which is in worse 
agreement with experiment than the three generation model. Thus, we conclude that independent of other exper- 
imental constraints, simply obtaining the observed CKM matrix structure for the three lightest generations in our 
model predicts that there is no fourth generation of quarks! 

We have also seen that the generic profile of the internal wave functions generates a rank three Yukawa matrix. In 
particular, this implies that the up quark is not massless, in accord with results from lattice gauge theory simulations 
of QCD. In particular, this implies that the strong CP problem is real. Interestingly, supersymmetry breaking in 
F-theory GUT models automatically contain an axion with a phenomenologically viable value for the decay constant 

0- . . 

We find it quite remarkable that the most minimal geometric structures in this framework are in fact sufficient 
for the purposes of achieving phenomenologically viable hierarchies in the Yukawa couplings. Turning the analysis 



9 



around, the beautiful match obtained in minimal F-theory GUT scenarios can be viewed as an important constraint 
on the class of compactifications, and in particular the requisite geometries which can reproduce detailed features of 
the MSSM. 

Finally, one can extend this analysis to the case of neutrino masses and mixing. Whereas the quarks and charged 
leptons localize on matter curves in 5'jin some F-theory GUT models the right-handed neutrinos iV/j localize on matter 
curves which only touch at a point [2] ■ We have argued that the Yukawa matrices for the quarks and charged leptons 
can exhibit similar hierarchies provided the two points of intersection p^p and Pdown are sufficiently close together. 
Assuming appropriate wave function overlaps so that the superpotential terms LNnHu and MmajNuNn are present 
will generate viable masses for the neutrinos via the seesaw mechanism. Note, moreover, that this type of geometry 
can also easily accommodate large mixing angles provided the interaction point leading to LEH^ on the L-matter 
curve is not finely tuned to be near the interaction point leading to LNjiH^. We are currently working out the details 
of this scenario. 

Acknowledgements The work of the authors is supported in part by NSF grant PHY-0244821. 

VI. APPENDIX: WAVE FUNCTION DISTORTION FROM FLUXES 

In this Appendix we provide additional details on the profile of wave functions in the presence of a background 
gauge field strength. We refer the interested reader to [2] for a more detailed analysis of wave function distortion 
effects due to gauge field strengths in the direction parallel to the matter curve. We present this example to illustrate 
how wave function distortion works in general. Indeed, it turns out that though gauge field fluxes alone do not distort 
the profile of the Yukawas, these contributions can mix with higher form potentials, leading to an appropriate Yukawa 
distortion (26i]. 

We first consider the case where this gauge field strength is constant, and then explain how the profile of the wave 
function changes in the presence of variations in this field. To be completely explicit, we consider the case of a bulk Eq 
theory with 27's localized along three matter curves of rank one enhancement to E-^ which form a triple intersection 
at a point where Eq enhances to E^- In the parent Eg theory, the background values of the scalar (f> and internal 
gauge fields satisfy the BPS equations of motion [li|: 

dAC^ = 0, = 0, F("'2) = 0, A F^^'^^ = 0. (53) 

The Eq theory is defined by Higgsing the parent Eg, theory in the direction: {(j>) — {(j>i)Ti + (^2)22 where the Ti 
denote two Cartan generators in the SU{3) factor of Eq x SU{3) C Eg. The loci of partial enhancement to E^ are 
then given by {(j>i) = 0, {(l>2) = and {(pi) + {(j>2) = 0. In addition, we also consider a background gauge field strength 
taking values in the Cartan of SU{3). 

Adopting a local system of coordinates zi, Z2 such that = {(f)i), the modes trapped along the curve zi = curve 
satisfy equations ([7]) and 

v2{di+Ai)^Pj+vi{d2 + A2)i^+z^Xi2 = (54) 
{dj + A-)xi2 + zi^i^^0 (55) 
{d^ + A^)xi2 + zii>2^0 (56) 

where in the above, we have taken a canonical presentation for the Kahler form. Similar equations hold for the 
other massless modes of the theory. When the background field strength is constant, the solutions to this system of 
equations depend to leading order on the fluxes through a term of the form cxp{A4^^ Zizj) , where A4^^ depends on the 
background field strength configuration. 

As an example, consider the special case where only F^j and are non-zero. A gauge field configuration which 
reproduces this field strength is: 

A ~ —F^jzi ■ dzi + F22Z2 ■ dz2. (57) 
Solving equations yields V'j — ^^d- 

ipj a X12 oc a (Z2) • exp (-zizT/ V^) ' exp (--^22 ' + (-^iT ' ^i^/2) + 0{y^)) (58) 

where a {Z2) denotes a holomorphic function of Z2. Adiabatically including the position dependence in the F^j generates 
a series expansion in the z's. 



Note added in revised version: 
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In the context of minimal F-theory GUTs, the 5 x 10 x 10 interaction term originates from a local enhancement 
from SU{5) to Eq. Here, the 5 localizes on a curve where SU{5) enhances to SU{6), and the lO's localize on curves 
where SU{5) enhances to 5'O(10). A subtle point in this construction is that when a full GUT multiplet localizes on 
the SO{lQ) curve, achieving a rank one matrix requires the lO's to localize on the same matter curve This appears 
to suggest the condition that the curve on which the lO's localize should self-intersect, or "pinch" at the point of 
enhancement. However, as pointed out in 25], the wave functions of the Q fields will now have two components, ^'g^ 
and , corresponding to the profile of the Q wave functions on the two local pieces of the pinched curve. Combined 
with the analogous contribution from the U wave functions, this would then lead to a sum of two Yukawa couplings. 
Because ^'g^ and ^q_ a priori have different coordinate dependence (as they localize on different pieces of the same 
pinched curve), this will generically lead to a rank two matrix, unlike the rank one case required in the analysis of 
this paper. A simple remedy to this issue is to impose a Z2 symmetry which interchanges these two components of 
the pinched curve. Alternatively, we can consider geometries quotiented by this same Z2 symmetry. In the quotiented 
geometry, the two enhancements from SU{5) to 5*0(10) are identified, with the Eq point invariant under this group 
action. More generally, one can consider geometries where a Z2 action interchanges two smooth curves. 

Although the existence of such Z2 quotients may appear ad hoc, this is in fact the generic situation in compactifi- 
cations of F-theory! As noted in 1], there could be monodromies acting on the seven-branes of F-theory. In fact such 
monodromies were crucial for understanding how the non-simply laced groups arise from F-theory compactifications 
[^ij. However, such brane monodromies were not used in the constructions of F-theory GUTs presented in fT*]. As 
pointed out in (25| , exactly these generic monodromies will make the Z2 quotient just described above automatic! 

The essential point is that in the breaking pattern Eq D SU{5) x SU{2) x C/(1), it is a rather special condition on 
the nature of the singularity to require two distinct 5*0(10) matter curves. Indeed, typically there will be a branch 
cut in the configuration, so that the two seemingly distinct components where a 10 localizes will interchange under 
monodromy around the point of enhancement. This corresponds to a single smooth curve with a branch cut locus 
emanating out from the point of enhancement. 

The relevance of such geometries for phenomenology was clarified in 25] where these configurations were analyzed 
as deformations of an element of the Cartan of Eq such that the Z2 Weyl group of the SU (2) factor in the above 
breaking pattern permutes the two lO's. This follows from the fact that the lO's transform as a doublet of SU{2). 
This is the same geometric Z2 quotient discussed earlier, and shows that rather than being special, it is in fact the 
more generic case from the perspective of geometry. 

It is important to note, however, that just as in any orbifold theory, working in terms of quantities invariant under 
the group action in the covering theory clearly suffices in all computations. Thus, it is always enough to compute 
all relevant wave function overlaps in the covering theory. Quotienting by the appropriate Weyl group (in the Eq 
breaking to 5C/(5) case, the Weyl group of 5C/(2)), it follows that the computation of Yukawas in the covering theory 
fully determines the Yukawas in the quotient theory, and does not modify the estimates presented in this paper. 

Finally, it has also recently been found in [2^] that background gauge field fluxes alone do not distort the rank of 
the Yukawa matrix. Nevertheless, the general philosophy of this paper that background fluxes of the compactification 
can alter the profile of the matter field wave functions, and also the structure of the Yukawa matrices was recently 
confirmed in [26j by considering a non-commutative deformation of the F-terms of the seven-brane superpotential 
induced by background -ff(i,2) = Hfj + tubH^s fluxes of the compactification. Such contributions correspond to 
higher dimension operators of the type briefly alluded to in equations ([7|) and ([8]) . In some simplified examples it was 
found in 26] that much of the hierarchical structure of the quark and charged lepton masses, and the CKM matrix 
are recovered. 
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